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We revisit the Kuramoto model to explore the finite-size scaling (FSS) of the order parameter 
and its dynamic fluctuations near the onset of the synchronization transition, paying particular 
attention to effects induced by the randomness of the intrinsic frequencies of oscillators. For a 
population of size N, we study two ways of sampling the intrinsic frequencies according to the same 
given unimodal distribution g{uf). In the 'random' case, frequencies are generated independently in 
accordance with g(u>), which gives rise to oscillator number fluctuation within any given frequency 
interval. In the 1 regular' case, the N frequencies are generated in a deterministic manner that 
minimizes the oscillator number fluctuations, leading to quasi-uniformly spaced frequencies in the 
population. We find that the two samplings yield substantially different finite-size properties with 
clearly distinct scaling exponents. Moreover, the hyperscaling relation between the order parameter 
and its fluctuations is valid in the regular case, but is violated in the random case. In this last case, 
a self-consistent mean-field theory that completely ignores dynamic fluctuations correctly predicts 
the FSS exponent of the order parameter but not its critical amplitude. 

PACS numbers: 89.75.-k, 05.45.Xt, 05.45.-a 


I. INTRODUCTION 

We revisit the Kuramoto model [1J to investigate the 
phase order parameter and its dynamic fluctuations near 
the onset of the synchronization transition. In partic¬ 
ular, we focus on their critical finite-size scaling (FSS) 
behavior. Even though the Kuramoto model and its ex¬ 
tended versions have been widely studied mm , some is¬ 
sues persist, in particular about the above points. In the 
few existing studies mm, the critical properties near the 
transition point have only been partially resolved, which 
prompted the present work. 

FSS is particularly important for data analysis, since 
real systems are always finite, and particularly so in many 
experimental systems, biological or not [Bj. Recently, the 
FSS of the synchronization order parameter in the Ku¬ 
ramoto model and several extended versions on complex 
networks has been investigated. An unusual FSS expo¬ 
nent v = | has been found HHg. This value was obtained 
in the classic case, hereafter referred to as the random 
distribution , where the individual frequencies of oscilla¬ 
tors are sampled independently according to a given uni¬ 
modal distribution g(to). In such a situation, oscillator 
frequencies can be arbitrarily close to each other even 
in a finite population. Interestingly, when the oscillator 
density fluctuation along the frequency axis is suppressed 
in a population with regular distribution , a very different 
FSS exponent v = | is obtained. Here we present de¬ 
tails of our numerical study that led to this previously 
reported value nfflnn- 


Dynamic fluctuations of the order parameter in the 
Kuramoto model is another long-standing problem [5] 
HU- In the perturbation theory developed by Daido 
it is claimed that they scale differently when approaching 
the synchronization threshold from above and from be¬ 
low, with the scaling exponent 7 = f in the supercritical 
region and 7' = 1 in the subcritical region. This result 
is quite surprising by itself, in the sense that conven¬ 
tional scaling theory predicts identical scaling behavior 
on either side of the transition. Note that Daido used the 
regular frequency distribution in numerical investigations 
as well as in his perturbation theory [4]. 

In this paper, we report extensive numerical results on 
the FSS of the order parameter and its dynamic fluctu¬ 
ation in the random and regular realizations of the in¬ 
trinsic frequency distribution. The FSS exponent in the 
regular case is estimated numerically as v ~ |, half the 
value of its random counterpart. In contrast to Daido’s 
predictions, we find that dynamic fluctuations scale with 
the same exponent 7 = 7' on both sides of the transi¬ 
tion. Interestingly, our data supports 7 ~ 7' ~ 1 for 
the random frequency distribution and 7 ~ 7' ~ | for 
the regular case. Given the order parameter exponent 
/3 = I, the hyperscaling relation 7 = S' — 2/3 is thus vio¬ 
lated in the random case, but obeyed in the regular case. 
We finally briefly discuss the implications of our obser¬ 
vations on the ordering process and dynamic response in 
the Kuramoto model. 

The paper is organized as follows. In Sec. II, we intro¬ 
duce the Kuramoto model and present the two types of 
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frequency realizations that we consider in this paper. In 
Sec. Ill, we investigate the FSS of the order parameter 
and its dynamic fluctuations with the random frequency 
distribution and discuss the validity of the hyperscaling 
relation in the model. In Sec. IV, we repeat the inves¬ 
tigation and analysis with the regular distribution. We 
summarize and discuss the results in Sec. V. Appendix 
A presents the mean-field approximation to the FSS in 
the case of the random frequency distribution. 


turns out that this quenched noise and the nonlinear dy¬ 
namic fluctuations are intimately related to each other in 
the sense that, when the former is removed, a very dif¬ 
ferent set of FSS exponents are obtained. The two types 
oscillator frequency distributions are discussed in Sec. Ill 
and IV, respectively. 


III. RANDOM FREQUENCY DISTRIBUTION 


II. THE KURAMOTO MODEL 


The Kuramoto model jl] consists in a finite population 
of N globally-coupled phase oscillators governed by 

,, K n 

j = 1, • • • ) AT. (2.1) 

fc=l 

Here <f>j is the phase of the j th oscillator, and uij its in¬ 
trinsic frequency. A given realization of the Kuramoto 
model corresponds to a particular choice of the N in¬ 
trinsic frequencies {tOj}. In this paper, we shall as¬ 
sume that the number density of these frequencies on 
the real axis follows, on average, the normal distribution 
g(w) = (l/v / 27r)e _ “ / 2 . The coupling strength K is pos¬ 
itive (K > 0). 

Phase synchronization is conveniently described by the 
complex order parameter defined as [T] 

1 N 

Z(t) = Ae ie = -^e^', (2.2) 

JV 3 =1 


where a nonzero (positive) value of A in the asymptotic 
N —► oo limit implies the emergence of phase synchro¬ 
nization, and 9 is the phase of the global order parameter. 
With this definition, Eq. (2.1) is rewritten: 


In this section, we consider the realizations (samples) 
of frequency sets {ox,} by randomly and independently 
drawing frequencies from the Gaussian distribution g(u>). 
We measure the phase order parameter and its dynamic 
fluctuation for each sample and average over many sam¬ 
ples. In practice, we numerically integrate Eq. (2.3), 
for system sizes up to N = 12 800, using the Heun’s 
method m with a discrete time step St = 0.01, up to 
total 2 x 10 5 time steps. We start with random initial 
values of {(p-j (0)} and average the data over time after 
the first 10 5 steps to avoid any transient behavior and 
also average over 10 2 — 10 3 different samples of {uij} and 
{<^•(0)}. At K = K c , we perform numerical integration 
for larger systems up to N = 204 800 with twice more 
realizations. 


A. Finite-size scaling of the order parameter 

In analogy with equilibrium critical phenomena, we 
expect the order parameter in the critical region to satisfy 
the FSS [14] . 


A(K,N) = N~ 0/p f(eN 1/D ) , (3.1) 

where e = K — K c and the scaling function f{x) is 
a monotonically increasing function with the following 
asymptotic properties, 


<f>j = uij — KAsm((/)j — 6 ), (2-3) 


where both A and 8 are dynamic variables defined by 
Eq. (2.2). In the limit N —> oo, the model exhibits a 
continuous symmetry-breaking transition at K = K c = 
2/[tt£K0)] = \/8jn PQ. In this limit, the phase order 
parameter is strictly zero on the subcritical side (K < 
K c ), while on the supercritical side (K > K c ), 


A ~ (K — K C )P with 



(2.4) 


This result holds for any unimodal and symmetric g(co). 

For a finite population of oscillators, A (f) exhibits both 
dynamic and sample-to-sample fluctuations whose scal¬ 
ing with N is the main focus of this work. The sample- 
to-sample fluctuations arise from the randomly drawn 
oscillator frequencies from the distribution g(u). Conse¬ 
quently, the spacing between neighboring frequencies on 
the real axis is Poissonian in any given realization. It 


/(a0 


or, x » 1; 

(—a;)^ - "/ 2 , x <C —1. 


(3.2) 


The exponent v is known as the FSS exponent. For large 
N , this scaling ansatz reproduces Eq. (2.4) on the su¬ 
percritical side, and yields A ~ N~P/ U at the transition 
and A ~ N~z in the subcritical regime, respectively. 

In the usual mean-field treatment of the Kuramoto 
model, the entrained state is established by solving 
Eq. ( |2.3| ) at a constant A and 9 (dynamic fluctuations 
are ignored): Oscillators with \ujj\ < KA get into an 
entrained (locked) state, reaching a fixed angle. On the 
other hand, the oscillators with \uij\ > KA perform pe¬ 
riodic motion with modified frequencies. With this, the 
self-consistent equation for A from Eq. (2.2), can be writ¬ 
ten as 


A = 



dujg{u>) \/\ — ( uj/KA) 2 


(3.3) 
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in the N —> oo limit. Eq. ( |3.3| has a nontrivial solution 
(A > 0) for K > K c = 2/ [ 775 ( 0 )]. 

For a hnite-size system, the right-hand side of Eq. (3.3) 
is replaced with 

A = T(KA) = 1 v/l-K/W, (3-4) 

which contains a sample-dependent correction dT = 
— [SE] oc y/Ng/N [8j. Here N s is the number of (en¬ 
trained) oscillators in the frequency interval (-K A, A'A) 
and [■] denotes sample average. As was shown previously 
in Refs. mm and explained in detail in the Appendix 
A, for K close to K c , Eq. (3.4) admits a scaling solution 
of the form <E3> with the exponents 



-1 -0.5 0 0.5 1 

e 


/3 = and v = \ ■ (3-5) 

The scaling function f(x), however, has significant 
sample-to-sample fluctuations. 

The exponent value j3 = \ is typical of various equilib¬ 
rium mean-field models described by the so-called the¬ 
ory. However, the FSS exponent value of v = | is quite 
unusual and implies much stronger fluctuations than that 
in ordinary models with the conventional value v = 2. 
This remarkable result has been confirmed by extensive 
numerical simulations BE! and the straightforward ex¬ 
tension to the Kuramoto model on various sparse random 
networks in ©muni. 

B. Dynamic fluctuations 


FIG. 1: (Color online) Dynamic fluctuation \ against t = 
K — K c for various system size N, with K c = y/WpH. The 
data have been obtained from the average over 100 — 1000 
samples with different sets of {o-'y} and initial phases {</b(0)}. 
Statistical errors are represented by symbol sizes. 

group (RG) sense, and hence 7 = 7 '. The corresponding 
FSS takes the form 

X(e, N) = AP'/ p /i(dV 1 '' p ) , (3.8) 

irrespective of the side of the transition, and the scaling 
function h(x) has the limiting behavior 

{ const, x = 0 ; 

a; -7 , i>1; (3.9) 

(-a:) -7 , x < — 1 . 


Now, we consider the dynamic (temporal) fluctuations 
of the order parameter by measuring the quantity m 

X(K, N)=N [((A - (A)) 2 )] = N [(A 2 ) - (A) 2 ] (3.6) 

where (•) denotes time average in the steady state of a 
given sample and [•] sample average, respectively. 

For K < AT C , we expect from simple power counting 
that x ~ 0(1), because A ~ N~ s. For K > K c , the or¬ 
der parameter becomes 0 ( 1 ), but these 0 ( 1 ) terms can¬ 
cel out each other exactly and the next order terms pro¬ 
portional to N~ 2 contribute. Again, \ becomes 0(1). 
However, at K = I\ c , strong critical FSS occurs since 
A ~ N~P/ lJ = N~s, and one can expect that x will 
diverge in the N —f 00 limit. Summing up, the critical 
behavior of x near the transition may be expressed, in 
the N —> 00 limit, as 


The two exponents 7 and 7 ' describe the divergence of 
X in the supercritical (e > 0) and subcritical (e < 0) 
regions, respectively. 

For homogeneous systems, scaling in the critical region 
is controlled by a single fixed point in the renormalization 


In this context, the Kuramoto model is somewhat spe¬ 
cial in that each oscillator has its own intrinsic frequency 
and the effect of the mean-field is different on different 
oscillators. Therefore the usual RG argument with only 
a few relevant variables may not apply and one should 
entertain the possibility that 7 / 7 '. In such a situ¬ 
ation, the FSS exponent v should also be different on 
the two sides of the transition, so that the ratio rela¬ 
tion 7 / 1 / = 7 'Iv' be satisfied, in order to maintain the 
continuity of x at e = 0. Indeed, through perturbative 
calculation, Daido [1] obtained 7 = \ and 7 ' = 1. 

We have carried out extensive numerical simulations 
to obtain accurate estimates of 7 and 7 ' as well as other 
FSS exponents. Figure |T] shows numerical data for \ 
versus e for various system sizes N. As expected, finite- 
size effects become huge near the transition point. We 
observe two important features of the finite-size effects: 
First, the peak position at K = K max is always on the 
subcritical side (e < 0 ) and approaches the bulk critical 
point (e = 0) as N increases. Second, the peak height 
Xmax grows rapidly with N. By analyzing the data for 
the peak position and height, we estimate 7 ' and v' in 
the subcritical side, using the conventional FSS ansatz 

Xmax ~ iV 7 ' /P ', <5K max ~ N-V*' (3.10) 
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FIG. 2: (Color online) (a) Behavior of Xmax is shown as a func¬ 
tion of N in a log-log plot. The straight line obtained through 
the least-square fitting method displays Xmax ~ N 0AQ . (b) 
Behavior of SK max (= |A' max — K c \) is shown as a function of 
N in a log-log plot. The straight line obtained through the 
least-square fitting method displays <5-Kmax ~ IV -0 ' 33 . 



FIG. 4: (Color online) Numerical data for the inverse of the 
effective exponent 7 ^ are plotted against e for N = 12 800. 
Two dotted lines represent the fitting lines, indicating 7 ss 
0.98 and 7 ' « 1.01. 

with 5K max = |.Km ax - K c \. As shown in Fig. l(a,b), 
both quantities exhibit power-law scaling with the expo¬ 
nent values 

V/z/= 0.40(1), l/v' = 0.33(8) . (3.11) 

Even though our estimate of l/v' is rather weak due to 
huge statistical uncertainties in locating K max , it includes 
our analytic value of 1 /v = | within the error bar. Tak¬ 
ing a simple ratio of the above exponent values, we esti¬ 
mate Y = 1.2(3). If we take l/v' = |, then 7 ' = 1.00(3). 

We also measure the value of x at the bulk critical 
point K = K c . As K c is located to the right of the peak 
position K max (N) for all N in this case, the FSS behavior 
at K c may be regarded as a continuation of the scaling 
on the supercritical side. Again, the FSS ansatz leads to 



N 


FIG. 3: (Color online) Critical increase of \ at K c (= y ^S/n) 
is plotted as a function of A in a log-log plot. The slope of 
the straight line is 0.39. Statistical errors are represented by 
symbol sizes. 


Xc = x(e = 0 ,N)~N->' P . (3.12) 

Figure [3] shows that the numerical data agree almost per¬ 
fectly with the FSS ansatz with 

7 /v = 0.39(2) . (3.13) 


This result confirms the ratio relation 7 /F = 7 ' /v\ but 
does not provide estimates for 7 and v separately. 

To estimate 7 and 7 ' directly from off-critical data in 
Eq. (3.7), we examine the local slopes of \ against |e| in 
double logarithmic plots for very large N. The effective 
exponent 7 e ff is defined as 


7 eff(e) 


din x 
din |e| ’ 


(3.14) 


and the critical exponents are obtained from their 
asymptotic values as 7 = lim e _ >0 + 7eff(e) and 7' = 
lim e _> 0 - 7 eff(e)- Figure 0 shows the inverse of the ef¬ 
fective exponent 7 A 1 as a function of e. We find that 
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7 e jf approaches y ' -1 linearly as e goes to zero from be¬ 
low in the subcritical region, with its asymptotic value 
7 ' = 1.01(3). This is consistent with our previous numer¬ 
ical estimates and also with the theoretical prediction of 
7 ' = 1 by Daido [J]. Note that 7 ' = 1 has been also 
obtained in a systematic 1/A expansion [5], where the 
random frequency distribution is prerequisite. 

On the supercritical side, varies in a nonlinear 
fashion far from the transition, but approaches 7~ 4 lin¬ 
early near small e with the asymptotic value 7 = 0.98(5). 
This suggests strongly that 7 = 7 ' = 1, and definitely ex¬ 
cludes Daido’s value 7 = With 7 = 1 , our numerical 


estimate in Eq. (|3.13|) lead to u = v' = §. 


Summing up all these results, we conclude that our nu¬ 
merical data strongly support simple scaling exponents: 


V= ~ V ' = \ 


and 7 = 7 ' = 1 . (3.15) 


However, a proper analytic treatment for the dynamic 
fluctuations in the supercritical side is lacking, and fur¬ 
ther investigations are thus needed. One should note that 
the dynamic fluctuations of the order parameter are not 
necessarily proportional to the susceptibility (response of 
the order parameter to an infinitesimal external field) in 
general nonequilibrium steady states due to probable vio¬ 
lation of the fluctuation-dissipation theorem. It is almost 
trivial to derive the susceptibility exponents in the Ku- 
ramoto model, which turn out to be 7 SUS = q( us = 1 [H]. 
Yet this accordance with the dynamic fluctuation expo¬ 
nents in Eq. (3.15) may be merely coincidental. 


C. Hyperscaling relation 


In equilibrium critical phenomenon, the equation 

l = (3.16) 


with v = du, is usually referred to as a hyperscaling re¬ 
lation. Here d is the dimension of the system and v is 
the correlation length exponent. Equation (3.16) follows 
from a simple power counting of E q. ( 3.6) at the tran¬ 
sition. The left-hand side of Eq. (3.6) is \c ~ A 7 ^, 
while the right-hand side is proportional to _/V 1_2/3 /" if 
the leading order terms do not exactly cancel out each 
other. 

It is obvious that our result Eq. (3.15) does not sat¬ 
isfy this relation. Violating the hyperscaling relation is 
only possible when the leading order terms in [(A 2 )] and 
[(A) 2 ] cancel out exactly and the sub-leading order terms 
show up, yielding 7=1 instead of |. A plausible con¬ 
jecture for leading and sub-leading terms for the large A 
expansion of the two quantities is then 


[(A 2 )] = A" 2/3 / p (c + dA-i) , 

[(A) 2 ] = r 2/3/p (c' + elV5) , (3.17) 



N 



FIG. 5: (Color online) (a) Critical decay of [(A 2 )] and [(A) 2 ] 
at K c (= y/R/'x) are plotted against A in a log-log plot. These 
two data are so close to each other that it is not possible to 
discern two data in this plot. The straight line is the fit¬ 
ting line of 1.14A -0 ' 40 . (b) As [(A 2 )] and As [(A) 2 ] versus 

A _ s. These data are expected to show straight lines. The 
least-square fitting yields two straight lines with the same in¬ 
tercept. Lines between data are only guides to the eye. 


with c = d. Indeed, our numerical data show in Fig. [ 5 ] 
that ck, d = 1.140(2), d = 0.10(4), and e = -0.04(4). 
Then, the dynamic fluctuation \ at the transition scales 
as 


Xc = {d- e)A 1 - 2/3 / p -s ~ A 7 / c 


(3.18) 


yielding 


7 = z/ —2/3—- = 1. 


(3.19) 


The cancellation of the leading order (c = d) in 
Eq. ( |3T7| would be expected if dynamic fluctuations 
only produce subleading order c orre ctions to the solu¬ 
tion of the mean-field equation ( |3.4| at the transition. 
Interestingly, as we show in Appendix A, the mean- 
field approximation yields [(A 2 )] = [(A) 2 ] ~ 0.93A _ § 
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which is smaller in amplitude as compared to the numer¬ 
ically determined value c = 1.14. Therefore it appears 
that dynamic fluctuations renormalize the strength of the 
quenched noise <54', but do not change its scaling with N. 

It is interesting to compare the violation of hyperscal¬ 
ing in the Kuramoto model to known cases in equilib¬ 
rium phase transitions. One well-known such case is the 
incompatibility between the correlation length exponent 
v = | at the Gaussian fixed point and the mean-field 
exponents /3 = ^ and 7=1 above the critical dimension 
do = 4 E|. Violation of hyperscaling is also reported 
in the three-dimensional random-field Ising model whose 
ordering transition is controlled by a zero-temperature 
fixed point [T 8 U 2 I] . In both cases, the energy associated 
with the symmetry-breaking field is much stronger than 
the thermal energy on large length scales, hence the order 
parameter fluctuations, as measured by the exponent 7 , 
become less than what is predicted by the hyperscaling 
relation. In this sense, the weaker dynamic fluctuations 
of the order parameter in the Kuramoto model may be 
interpreted in a similar way, despite the fact that we are 
now in a nonequilibrium situation. 


IV. REGULAR FREQUENCY DISTRIBUTION 


In this section, we consider the ‘regular’ frequency set 
{ojj} with minimal disorder among frequencies. This 
set can be generated following a deterministic procedure 
given by 


j - 0.5 
N 




(4.1) 


with j = 1 and g(co) = (l/v / 27 r)e -aj2 / 2 . The 

generated frequencies are quasi-uniformly spaced in ac¬ 
cordance with g(uj) in the population. As this set is 
uniquely determined for each N, sample disorder may 
result from initial conditions only. Note that the regu¬ 
lar frequency distribution is considered in numerical at¬ 
tempts by Daido [4]. 

Similar to the case of the random distribution, we 
numerically integrate Eq. (2.31 for system sizes up to 
N = 12 800 with a discrete time step St = 0.01, up to to¬ 
tal 10 7 time steps. We start with random initial values of 
{^( 0 )} and average the data over the latter half of time 
steps and also average over 20 — 100 different samples 
of {c/(0)}. At K = K c , we perform numerical integra¬ 
tion for larger systems up to N = 25 600 with 10 — 100 
samples. 


A. Finite-size scaling of the order parameter 

We measure the order parameter and its square dur¬ 
ing numerical integrations and average over time in the 
steady state, and then average over samples with differ¬ 
ent initial values for the phase variables. These averages 



N 

FIG. 6 : (Color online) Critical decay of the order parameter 
[(A)] and [(A 2 )] at K c (= 1/8/71-) is plotted as a function 
of A in a log-log plot, which are represented by (red) open 
boxes and (blue) open circles, respectively. The slopes of two 
straight lines are —0.39 and —0.78, respectively. 



FIG. 7: (Color online) Scaling plot of the order parameter A 
for various system size N for the regular distribution. We use 
P/v = | and v = |. 


are denoted by [(A)] and [(A 2 )], respectively. Figure krt 
shows their size dependence at the transition (e = 0)7 
exhibiting the nice power laws 

[(A)] ~ at-°' 39(2) and [(A 2 )] ~ AT-°' 78(3) . (4.2) 


Following the FSS theory in Eqs. and ( |3.2| ), we es¬ 
timate the decay exponent as 


/3/P = 0.39(2) 


(4.3) 


Given the value /3 = \, we obtain P ss | . 

As a cons istency check, we investigate the FSS relation 
in Eq. (3.1) for e / 0. We plot AIV^ 1 ' versus eN 1 ^ for 
various system sizes N in Fig. [71 (For a moment, we drop 
the average brackets as A = [/A)] for convenience.) The 
data show a perfect collapse with the choice of P = |. We 
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FIG. 8 : (Color online) Dynamic fluctuation \ is plotted as 
a function of e(= K — I \ c ) for various system size N for the 
regular distribution, where K c = y/8/iv. The solid (red) line 
represents a fitting line for data with largest system size. 


FIG. 9: (Color online) Critical increase of \c at K c (= y/8/Tr) 
is plotted as a function of N in a log-log plot for the regular 
distribution. The slope of the straight line is 0.22. 


also considered the regular Lorentzian distribution given 
by g{ui) = 7cj/[7r(w 2 +7 w 2 )] with the half-width 7 ^, where 
the oscillator frequency is chosen as u>j = 7 ^ tan[j7r/iV — 
(N + l)7r/(2JV)] for j = We find a similar 

behavior again with u m |. 


B. Dynamic fluctuations 


In this subsection, we present numerical data for the 
dynamic fluctuations x i n Fig. [ 8 ] Compared to the ran¬ 
dom distribution case, \ is very much reduced in magni¬ 
tude and peaks only in a much narrower region around 
the transition. The small amplitude makes it difficult 
to determine the peak height and position accurately. 
However, with the exact information of K c = -y/ 8 / 7 r, it is 
numerically feasible to investigate the critical increase of 
Xc at the transition as a function of N. The data shown 
in Fig. [9] fit well to a power law 

X c ~ N^/ v with 7 /v = 0.22(2) . (4.4) 


Using v = |, we find 7 = 0.27(3), which is close to 
This value is consistent with Daido’s theoretical result of 

7= 3 I3- 

We also estimate 7 and 7 ' directly from off-critical data 
by exa min ing the effective exponents 7 e fi for large N. 
Figure 10 shows ( 47 e ff ) -1 versus e. Following a similar 
procedure as in the random case, we find 7 = 0.24(2) 
and 7 ' = 0.25(1). We also fit the off-critical data directly 
in Fig. [ 8 j which agree very well with \ ~ |e| —0-25 . All 
the numerical results strongly suggest scaling exponent 
values as 




and 7 = Y = - 
4 


(4.5) 


Our results are clearly different from Daido’s numerical 
suggestions of v ss 4 and v' « 2 and also his theoret¬ 



FIG. 10: (Color online) l/( 47 e ff) versus e for N = 6400. Two 
dotted lines represent the fitting lines, indicating 7 fts 0.24 
and x' ~ 0.25. 


ical prediction of 7 ' = 1 on the subcritical side, which 
cannot be accounted for by adopting a slightly different 
definition of \ by Daidopj)]. 


C. Hyperscaling relation 

The hyperscaling relation, 7 = v — 2/3, is seen to hold 
in the regular case. Furthermore, there is no detectable 
“disorder fluctuation” arising from different initial phase 
values, leading to [(A )] 2 « [(A) 2 ]. 

V. SUMMARY 

In this paper, we revisited the well-known Kuramoto 
model and investigated the FSS of the order parameter 
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and its dynamic fluctuations near the onset of synchro¬ 
nization transition. We find that critical scaling behav¬ 
ior crucially depends on the presence of disorder among 
intrinsic frequencies of oscillators. In the random case 
with disorder, it is shown analytically and numerically 
that the FSS exponent v = | and the dynamic fluctua¬ 
tion exponent 7 a; 1 , in both the supercritical and the 
subcritical region. In the regular case without disorder, 
we find v ss | and 7 « In both cases, the exponent 
v is different from its conventional value v = 2 in usual 
mean-field theory of homogeneous systems with global 
coupling. 

In the random case, disorder in the oscillator frequen¬ 
cies broadens the critical region and drives order param¬ 
eter fluctuations in a sample-dependent manner. The hy¬ 
perscaling relation that describes the interplay between 
static (i.e. averaged) and dynamic order parameter fluc¬ 
tuations is violated. Interestingly, a modified mean- 
field theory that completely ignores dynamic fluctuations 
yields correct critical FSS of the order parameter, al¬ 
though the amplitude is renormalized by dynamic fluc¬ 
tuations. This is one of the main findings in this work. 
In comparison, numerically determined exponents 7 and 
v in the regular case do satisfy the hyperscaling relation. 
Further work is required to explain the particular values 
observed in our study. 
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Appendix A: Finite-size scaling in the mean-field 
approximation 


In this appendix we consider in some detail the solution 


to the self-consistent mean-field equation (3.4). For this 


purpose, it is convenient to introduce a function 


ip(u) = 




1 0 , 


< i; 
> 1 . 


(Al) 


In terms of -0, we have, 


N 


*(KA) = -Y / ^j/KA). 


3 =1 


(A2) 


Since each term in the sum is a monotonically increasing 
function of z = ILA, independent of the sign of Uj, the 
sum is also a monotonically increasing function of z. In 
terms of z, Eq. (3.4) can be rewritten as, 


(A3) 


When I\ is sufficiently small, only the trivial solution 
z = 0 is obtained. The nontrivial solution is obtained 
only for K > K c where K c depends on the specific choice 
of the frequencies {ujj}. 

We now co nsid er the statistical properties of the solu¬ 
tion to Eq. (A3) when the frequencies {iOj} are drawn 
independently from a given distribution g{u>). In this 
case, ijjj = ip(u}j/z) are also independent random vari¬ 
ables whose distribution can be obtained from g(ui). The 
quantity 4/ as defined by Eq. (A2) is simply the average 
of the random variables tpj, j = 1,..., N. Its mean value 
is given by [ip\. Let (5'F = 'F — [ip]. Its variance is given 
by 


[(<5T) 2 ] = 


m - m 2 

N 


(A4) 


The variance of ip can be easily computed from g(u>): 

[ 1 P 2 } - [ 1 P ] 2 = J dug(u)(l- (^) 2 ) 


-(l^sM^i-^) 2 )' 

g(o)z- ( 


(A5) 


Invoking the central limit theorem, we conclude that the 
distribution of *F tends to a Gaussian at large N, with 
mean and variance as stated above. 

For small z, 


= J dug{u})^j 1 - (^) 2 

= |-7(0) 2 +^ff"(0)z 3 + ... 


Equation (A3) can now be written in the form 


K z = K„ z — cz 6 


(A 6 ) 


(A7) 


where K c 1 = ng{Q)/2, c = — 7rg"(0)/16 > 0, and d = 
(4g(0)/3)5. Here 7 = d’F/fd’F 2 ]® is a Gaussian random 
variable with zero mean and unit variance. 

Sol ution to Eq. 1 A7h can be cast in the scaling 


form (3.1) with /? = 2 an d v = Specifically, let 
z = K C N~^ /, Eq. (A7) becomes 


xf-cK%f 3 + = 0, 


(A 8 ) 


where x = N b (K — K c ) /K is the scaled distance to the 
transition point. 

We now consider solution to Eq. (A 8 ) at K = K c or 
x = 0. Simple algebra yields 


/( 0 ) = 


Agl , 77 > 0; 
0, 77 < 0. 


K 1 z = \F(z). 


(A9) 
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Here A = (8/37 r)sc sA' c 5 . The nth moment of /(0) at 
the transition is given by (n > 0), 


irm 



(A10) 


In the special case g(u) = e “ 2 / 2 /VTn, A = (87r/3)e. 


The first four moments are then given by 


[/( 0 )] 
t/ 2 (0)] 
[/ 3 ( 0 )] 
[/ 4 ( 0 )] 



0.643454.., 

0.930853.., 

1.45621.. , 

2.41398.. . 
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